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Abstract 

We investigate the coherence properties of pairs of counter-propagating atomic clouds, produced 
in superradiant Rayleigh scattering off atomic condensates. It is shown that these clouds exhibit 
long-range spatial coherence and strong nonclassical density cross-correlations, which make this 
scheme a promising candidate for the production of highly directional nonclassically correlated 
atomic pulses. 
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I. INTRODUCTION 



Atom-atom correlations, and the engineering of "nonclassical" states of atoms that exhibit 
quantum correlations, are currently attracting intense theoretical and experimental inter- 
est in the framework of ultracold quantum gases. In most of these studies, Bose-Einstein 
condensation plays a central role, since correlated atomic beams can be generated by ma- 
nipulating appropriately preformed condensates Directionality of the produced atom 
pairs, is of vital importance for many potential applications (e.g., subshot noise precision 
measurements, tests of quantum nonlocality, and presumably quantum information process- 
ing) but, as discussed in js], it is difficult to be achieved within existing schemes without 
seeding (e.g., see 



Superradiant Rayleigh scattering from an elongated atomic condensate 



9|, has the po- 



tential to produce highly directional counter-propagating matter waves, which exhibit cor- 



relations as shown in 



10|. Moreover, in contrast to other proposals, the two matter waves 



have well defined spatial profiles, and one can tune their macroscopic populations. So far, 
however, the usefulness of the scheme remains debatable, since there has been no theoret- 
ical quantitative analysis of the coherence properties of the matter waves, and the type of 
correlations involved. 

Most of these issues are addressed in this article, within a model that treats the scattered 
photons and the matter waves quantum mechanically. Our model is expected to describe 
accurately the essential aspects of the process. Moreover, it includes spatial propagation 



effects, which are crucial 
and were not included in 



br the thorough understanding of condensate superradiance 
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111, 



12|. 



It is shown that the counter-propagating atomic clouds exhibit long-range spatial coher- 
ence. Their density cross-correlations are enhanced due to the formation of correlated atom 
bunches, which are also responsible for violation of the classical Cauchy-Schwarz inequal- 
ity, relative number squeezing and entanglement. Unlike earlier schemes js-?!, we study 
correlations in the context of mixing two optical- and two matter waves, hence the name 
"Four Heterowave Mixing". The present results shed light on the physics of such a mixing 
process, and they determine the operational regimes for possible applications and future 
experiments. 
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II. MODEL 



The system involves an elongated condensate of length L, oriented along the z axis and 
consisting of atoms jsl. It is exposed to a linearly polarized pump laser pulse E;(x, t) = 
Ei{t)G^y\&'^^^^~'^^*^ + c.c.]/2, with uj[ = cki, traveling in the x direction. The laser is assumed 
to be far off-resonant from any atomic transition, inducing thus Rayleigh scattering. 

Due to the coherent nature of the condensate, successive Rayleigh scattering events are 
strongly correlated, leading to collective superradiant behaviour |9l-|l2[|. Moreover, as a 
result of the cigar shape of the condensate, the gain is largest when the scattered photons 
leave the condensate along its long axis, in the so called endfire modes with momenta 
k ^ ^h^z and frequency u ^ ui. As a consequence, the recoiling atoms have well-defined 
momenta and appear in distinct atomic side modes. In the side mode (n, m), atoms possess 
momentum q ^ hki {ne^ + me^ ) , while the corresponding frequency is approximately given 
by uJn,m ~ {n"^ + m'^)ur, where Ur = hkf /2M is the recoil frequency, and M the atomic mass. 
In this notation, the "side mode" (0,0) describes the condensate at rest, while (1,±1) and 
(— 1, ±1) are the first-order forward and backward atomic side modes, respectively. 

We are here concerned with the early stage of the process, where only first-order atomic 
side modes become significantly populated. The condensate remains practically undepleted 
and can be treated as a time-independent classical field $(x) = ip{z)ipj_{x, y), where ^p{z) and 
f±{x,y) are the longitudinal and transverse wave functions. For short times, the coupling 
between the counter propagating optical end-fire modes can be neglected. Taking advantage 
of the symmetry of the system with respect to the x-axis, we focus on the dynamics of the 
single endfire mode (taken to be monochromatic with k = +kiez), and the first-order side 
modes coupled by it [i.e., (1, —1) = (+) and (—1, 1) = (— )]. 

Under the slowly- varying-envelope approximation (SVEA), we decompose the opera- 
tors for the matter wave and the positive- frequency electric fields as \l/(x, t) = $(x) -|- 
^±J2j^ji^,t)(i~'^'^'^~^^'"''~^'^\ and E(+)(x,t) = E{+^(x,t) + ^J^u±e{z,t)eye-'^'^^~^'\ with 
u±{x,y) the transverse endfire-mode profile, and j G {+, — }. The operators ipj{z,t), e{z,t) 
are annihilation operators for side- and endfire-mode fields with transverse behavior fixed 
by f± and u± (see appendix). By construction, these operators obey standard commuta- 
tion relations with, however, Dirac delta functions having finite width of about 1/k (see 
appendix). Inserting the SVEA expansions into the Maxwell- Schrodinger equations for the 
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coupled matter- wave and electric fields 13| , one can derive equations of motion for the oper- 
ators 'ijjj{z,t) and e{z,t) Rescaling to dimensionless time r = 2ujrt and length ^ = kiz, 
we obtain 



i^^i(e,r) = -«:e(e,r)^*(0, 
i^^Ue, r) = Kei^, tMO + r), 



dr 



(la) 
(lb) 

(Ic) 



with = ?/'_e~^''^, X = cki/{2uJr) and k, = g\/kiL/{2ujj) [15|. The atom-photon coupling is 
given hy g = \d\'^£i/ {2h^5) ^Jhw / {2eQL) J dxdy ip\u±, where d is the atomic dipole moment, 
and 6 the detuning of the laser from the nearest atomic transition. By discarding back- 
wards recoiling atomic modes in Eqs. ([T]), the remaining equations describe conventional 
superradiance. 

Solutions to the system of Eqs. ([1]) can be expressed as integrals involving the operators 
evaluated at the boundary of their domain — i.e. at ^ = and r > or vice versa at 
^ > and r = — by applying Laplace transform techniques. Due to the large value of x, 
retardation effects can be neglected, and the solutions for the matter wave operators read 



V^Ue,r) = V^Ue,0) 

+ inip*iO I rfr'e(0,r')Fi,o(75,o,r-r 



(2a) 



+ r¥^*(0|y^ ^eV(0^i(e',o)F2,o(7a',^) 

'%eV*(e')^-(e',0)Fi,i(75,e,r)|, 

- iKifi^) I dr'e(0,r')Fo,i(75,o,r-r') 

- r^(o{^ rfev(e>l(r,o)Fi,i(75,e,r) 

+ ^^rfeV*(e>-(r,0)Fo,2(7c,e>^)}- (2b) 
We have introduced V = n'^/x and the functions F^^^{a,f3) = 
C~\p \^e°'^^e~°'^^'^'^^^p^^{p + 2i)~'^^ , with C^^ denoting the inverse Laplace transform. 



while 7^^^' = r J^, d(\ip{()\'^ . The functions F^^^ can be expressed exphcitly as combinations 
of Bessel functions and integrals thereof (see appendix). 

Experimental observations of superradiance from condensates distinguish between the so- 
called weak- and strong-pulse regimes jol- Typically, the former (WP) regime is characterized 
by coupling constants g ~ 10^s~^ and TN ^ 1, while for the latter (SP) regime g ~ 10®s~^ 
and FA^ ^ 1. Throughout our simulations, we chose FA^ = 1 and FA^ = 100 for the two 
regimes. The condensate at rest was taken to be Thomas-Fermi distributed, i.e. f{z) = 
^/'oiLz — z'^)Q[z) / , with L = 130/im and &{■) the Heaviside step function; containing 
iV = 10^ ^''Rb atoms. The incoming laser pulse was chosen to have a rectangular profile and 
ki = 8.05 X lO'^m-^ 



III. COHERENCE AND CORRELATIONS 

The first-order correlation function of the side mode (j) is G*^j^(^i, ^2; t) = 
Wj{^i,T)ipj{C,2,'T~)) ■ It characterizes the coherence properties of the side mode, and quan- 
tifies the characteristic coherence length over which phase correlations exist. Using Eqs. 
(E]), the commutation relations for the fields [l6|, and applying the boundary values 
(e^e(,^ = 0,r)) = {'4^j'4'j{^,T = 0)) = one can easily obtain analytic expressions for 
for any time r. In practice, however, it might not be feasible to resolve the 
exact atomic positions .^1, and .^2 in the two sidemodes. For direct comparison to possible 



experiments, we consider the volume averaged degree of coherence, defined as 17|] 

where Ar,y(ei, A^; r) = {nj{^,, r)) (%,(ei+Ae, r)), with r)) = r)4-(e, r)), = 

,^2 - 6 and A = kiL. 

As depicted in Fig. [H g^^\AC,; r) attains its maximum for A^ = 0, where gjf{0; r) = 1. 
Let the coherence length of the side mode (j) be the spatial separation Xj for which, 
gj^\Xj;T) ~ 1/e 18|. The matter-waves associated with the two side modes can be con- 
sidered as first-order-coherent for any two points separated by |A,^| < Xj. According to 
Fig. [U the long-range order of the condensate is transferred to the atoms in the side modes 
(±). More precisely, for relatively short times, the coherence length of the side mode 
is Ai ~ A, while for the side mode (— ) it is somewhat smaller A2 ~ 0.6A. This is because 



FIG. 1: Degree of first-order coherence for the j'th side- mode at various times, in the WP (a) and 
SP (b) regimes. The function is symmetric with respect to = 0. 



backward recoiling atoms are created when incoherent Ught of the endfire mode is scattered 
by condensed atoms (as opposed to forward recoiUng atoms, which involve scattering of co- 
herent laser light), and thus the transfer of coherence is not as efficient as for the side mode 
(+). In the SP regime, the coherence length of the side mode (+) decreases as time goes on, 
approaching the one of the backward side mode, which remains practically constant. In the 
WP regime, however, the corresponding drop is not so prominent, and instead we observe a 
growth of the spatial coherence for the side mode (— ). 

Density correlations are described through the second-order correlation functions 
Gfji{ii,^2]i') = f>'j{^i,T)hji{^2,T) '■)■ This quantity reflects the probability of finding 
an atom of type (j') at a position ,^2, given that an atom of type (j) has been detected at 
position ^1. Using Eqs. ([2]), one can express G^y(^i,^2;T) in terms of first-order correla- 
tion functions. The corres pon ding volume averaged normalized correlation functions can be 
defined in analogy to g^^^ 17|, obtaining 



Jo d^iAljf{^i,A^;T) 



;;(2)/Ac.^^ _ 1 , Jo 



for pjj = Gfh^i,^2]T) and p+_ = (^'+(6, r)^/'_(^2, r)). 
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In view of Eq. ([3]), the function Gjj {^i,$,2',t) determines to a large extent the density 
autocorrelation function gj^\A^;T). The behavior of ^j^^(A^;r) throughout the evolution 
of the system resembles the behavior of the first-order coherence, albeit with somewhat 
different profiles, not shown here. In contrast to ^j^^(A,^;r), however, ^j^^(0;r) — )■ 2, as 
|A^| —7- 0. This is a manifestation of correlations between atomic densities at two different 
points separated by |A^| <^ A^. In other words, the recoiling atoms in the side mode (j) 
tend to appear bunched and thus, detecting an atom in position ^i, significantly increases 
the probability of detecting another atom close to it. The spatial extent of atom bunches is 
of the order of A-,-. For |A,^| > A-,-, gjj {0]t) — )■ 1, indicating that atomic densities become 
uncorrelated for larger separations. 

Density cross-correlations between the two side modes, are described through g_^__{AC^, r), 
which is plotted in Fig. [2l For any time r, there is a discontinuity at = ^2 due to one 
of the terms in p+-, which involves — C,2)- The physical reason for this discontinuity 
is the fact that an atom at ^1 in the sidemode (-I-) can be correlated to another atom at ^2 
in the sidemode (— ), only when the latter is created by scattering the endfire photon that 
was emitted by the former. This process (to be referred to as photon exchange hereafter) is 
only possible for ^2 > ^1 since endfire photons possess momenta +/cez. 

The enhanced correlations observed in Fig. [2] for |A,^| < A+, arise through the interplay 
between atom bunching, exhibited by the side mode (+), and photon exchange. A recoiling 
atom of type (— ), is not correlated only to the atom of type (+) that has emitted the photon, 
but rather to the entire bunch. Moreover, due to bosonic enhancement, the production of 
a backward-recoiling atom at a given position, stimulates the creation of additional atoms 
of the same type nearby. We have thus the production of correlated atomic bunches of type 
(+) and (— ). This is clear in Fig. [21 where the cross-correlation function is always larger 
than unity, and increases with decreasing |A^|, attaining its maximum value for A,^ — )■ 0^. 

Thus, measuring an atom in one mode significantly increases the probability of measuring 

(2) 

one in the counter-propagating mode. The spatial dependence g\_L{A^,T) at various times, 
bears analogies to the corresponding behavior of g^^\AC,, t). As a general remark, in the SP 

(2) 

regime, we find a much faster decay of (7^„(A^, r) compared to the WP regime. 

A further, qualitative analysis of correlations can be obtained by means of inequal- 
ities. It is known that "classical" fields satisfy the inequality V(^i,^2;t) < with 
V = - for all times r and pairs (^1,6) Q- 
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(a) 








(b) 





FIG. 2: Normalized density cross-correlation function at two different times, (a) WP regime: 
T = 1.025 (+), r = 4.025 (□); (b) SP regime: r = 0.03025 (+), r = 0.12075 (□). 
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FIG. 3: (Color Online) (a,b) Density plots of V at TNt ~ 3. (c) Time evolution of S (upper) and 
£ (lower). WP regime: (b), dashed and dot-dashed lines; SP regime: (a), solid and dotted lines. 

A quantum field, however, can violate this inequality if its P-representation attains negative 
values in a certain region of space and/or time. As depicted in Figs. |3t^a,b), in our system, 
photon exchange may lead to violation of this inequality only for .^2 ^ Such a violation 
can be also associated with squeezing or entanglement. The squeezing in the population 
difference between the two side modes is quantified by, S{t) = 1 + {R+)^^{: (Ai?__)^ :), 
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where AO = O - (6), R±{t) = N+{t) ± iV_(r) and iV,(r) = J^d^fiji^^r) As depicted 
in Fig. [3]^c), the system exhibits squeezing in the SP regime only where S{t) < 1. This 
never happens in the WP regime where (: (A_R„)^ :) ^ /J^ /J^ c^'Ci'^'C2|G'!,^^({i, ,^2; t)^. A 
necessary condition for separabihty of the bipartite system consisting of the mod es ( ±) is 



0, where £(r) = var(i+ + i_) + var(fi+-fi„) and Q = \{[Aj,Bj])\ (l9j. We 
have checked this inequahty for Aj^r) = Jq d^i'ipj + and Bj^t) = i d^{ipj "V^]); which 
are the quadratures of the side mode (j), with respect to its mean momentum. As depicted 
in Fig. [3]^c), the two side modes are indeed entangled in the SP regime for < FA^r < 4.2. 
When £^ > 0, the above criterion does not allow us to infer anything about the entanglement 
in the system. In view of these results, the SP regime seems to be more interesting for 
practical purposes than the WP regime. 

IV. CONCLUSIONS 

Superradiance from condensates is a promising technique for producing atomic clouds, 
with well defined momenta, spatial profiles, as well as tunable macroscopic populations. 
The counter-propagating matter waves typically produced in the strong-pulse regime, ex- 
hibit long-range coherence and nonclassical correlations, which can be explored for various 
applications. Our description does not take multi-modal superradiant emission into account 
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1C|, which can be strongly suppressed by reducing the aspect ratio of the condensate 



13|. The main predictions of our model are expected to be valid in more elaborate 
three-dimensional simulations, which can take into account also the propagation of atoms. 
In that case, however, the discontinuous gap in g\_{^^:, r) is expected to appear as a sharp 
but continuous increase. The experimental techniques that have been developed over re- 



and should be 



cent years allow for the direct measurement of correlation functions 
applicable to the verification of the present theoretical predictions. Finally, some of the 
present results might also apply to other systems which exhibit similarities to superradiance 
off condensates such as the collective atomic recoil laser 

The work was supported by the EC RTN EM ALL 
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2l|. 
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Appendix A: Operator Expansion 



Our operators are expanded as 



TT 



peAo 



(Al) 
(A2) 



peAo 



where j E { — 1,+1}, UJ±i 



2M 



and Aq is an interval around zero in momentum space 



which cannot be chosen larger than {—k/2, k/2) for ip+{^y 0); i^-i^', 0) 



to hold. There- 



fore the equal field commutator ipjiz, 0),ipj{z', 0) 
a distribution with width of the order 1/k. 



is not a Dirac delta function, but rather 



Appendix B: Inverse Laplace Transforms 



Fi,i{y,z 
F2,o{y,z 

^0,2 (y, 2 



-i2z' 



F,,o{y, z) = Jo {2^) Q{z) - e{z)^ "^-"'T? ^0 [2^/y{z - z')] J, (2 (Bla) 
Fo,i{y,z 



-i2z' 



Jo \Jz - z' 
e(^)^ dz'e-'^^'h \^^y{z - z')^ Jo {2^ 



h 



2v^ 



z — z' 



Jo[2^yi^lh) 
(Blc) 



-h{2^)Q{z)-Q{z) / dz'e 

y Jo 



I -2\z' 



Z — Z 



—h [2 ^y{z - z')] Ji [2 (Bid) 



e-''\ -Ji{2JlJ^)Qiz) + Qiz) / dz'e 

y Jo \ z- z' 



-h[2y/yiz-z')]M2^W^) 



with Jj and Jj the ith Bessel function of the first kind and the ith modified Bessel function 
respectively. 

Appendix C: Correlation Functions 

+5,-+M(ei,6) / rfr'Fi,o(7Ci,o,r-r')F*o(76,o,^-^'), 
Jo 

(CI) 
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where M(^i,^2) = ^f{Ci)f{^2), is Kronecker's delta, and S = min(^i, ^2)- 



G ■l(6,6;r) = (%(ei,r))(n,v(6,r)) + |a,y(6,e2;r)r, 



(C2a) 



where 



(6,e2;r) = G$.;.H6,e2;r), 



(C2b) 



and 



(V;+(ei,r)V;_(6,r)) 



M(ei,6)i^i,i(76,en^)0(6-ei) 



Jo 



(C2c) 
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